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ABSTRACT 

We present first results from the Southern Cosmology Survey, a new multiwavelength survey of 
the southern sky coordinated with the Atacama Cosmology Telescope (ACT), a recently commis- 
sioned ground-based mm-band Cosmic Microwave Background experiment. This article presents a 
full analysis of archival optical multi-band imaging data covering an 8 square degree region near right 
ascension 23 hours and declination -55 degrees, obtained by the Blanco 4-m telescope and Mosaic-II 
camera in late 2005. We describe the pipeline we have developed to process this large data volume, 
obtain accurate photometric redshifts, and detect optical clusters. Our cluster finding process uses 
the combination of a matched spatial filter, photometric redshift probability distributions and rich- 
ness estimation. We present photometric redshifts, richness estimates, luminosities, and masses for 
8 new optically-selected clusters with mass greater than 3 x 1O 14 M at redshifts out to 0.7. We 
also present estimates for the expected Sunyaev-Zel'dovich effect (SZE) signal from these clusters 
as specific predictions for upcoming observations by ACT, the South Pole Telescope and Atacama 
Pathfinder Experiment. 

Subject headings: cosmic microwave background — cosmology: observations — galaxies: distances 
and redshifts — galaxies: clusters: general — large-scale structure of universe — 
methods: data analysis 



1. INTRODUCTION 

The new generation of high-angular resolution Cos- 
mic Microwave Background (CMB) ground-based exper- 
iments repres ented by the the Atacama Cosmolo gy Tele- 
scope (ACT) ( |Kosowsky|2006l|Fowler et al.|2007D and the 
South Pole Telescope (5PT) ( |Ruhl et al.||2004| ) are cur- 
rently targeting their observations m a common area in 
the southern sky that will ultimately cover several hun- 
dreds to thousands of square degrees. These experiments 
will provide a blind survey of the oldest light in the Uni- 
verse at wavelengths of 1 — 2 mm and angular scales 
beyond the resolution limits of the WMAP and Planck 
satellites. At these arcminute angular scales, temper- 
ature fluctuations in the CMB are dominated by sec- 
ondary effects arising from the formation of large-scale 
structure in the universe. One of the strongest effects is 
the imprint left by galax y clusters though the Suny aev 
Zel'dovich effect (SZE) ( ISunyaev fc Zeldovich|[T980l ) in 
which CMB photons suffer inverse Compton scattering 
by the hot intracluster gas. ACT and SPT are designed 
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to detect the SZE, through its frequency-dependence: 
these experiments will measure temperature shifts of 
the CMB radiation corresponding to a decrement be- 
low and an increment above the "null" frequency around 
220 GHz. 

Much can be learned about the Universe from these 
surveys. First, accurate systematics-free maps will allow 
measurement of the primary power spectrum of tempera- 
ture fluctuations at all scales on which they are the dom- 
inant contribution. Second, these data sets will result 
in a complete census of massive clusters to arbitrarily 
large distances, limited only by a minimum cluster mass 
set largely by the instrumental sensi tivity and expected 
to be several 10 14 M (|Ruhl et al.1|2004| |Sehgal et al 



2007). Thanks to the relatively clean selection function 



as well as the redshift independence of the SZE, the clus- 
ter sample, especially the evolution of the number den- 
sity of clusters with redshift, will be quite sensitive to 
the growth of structure in the Univers e offering a po- 
tentia lly powerful probe of dark energy QCarlstrom et al. 
2002). Moreover, the SZE data in combination with op- 
tical, UV and X-ray observations can teach us a great 
deal about the detailed physics of cluster atmospheres 
and galaxy evolution in these dense environments. 

Significant observing time and effort has been devoted 
to the development of techniques and the detection of 
galaxy clusters using large-area optical catalogs and X- 
ray observations. Several projects have taken advantage 
of large-area CCD imaging and have developed auto- 
mated cluster detection sche mes to produce large cata- 
logs of clusters of galaxies (see iKoester et al. 1120071 iPost- 
man et al.|[T99 6, 2001 ; [Gal et al.|2000l | 2003l ( 2009 T|UCT 



ders &; Yee||2005| for example) which target the relative 
over- abundance of galaxies over a range of redshifts. Sim- 
ilarly, X-ray surveys such as the R OS A T All Sky Survey 
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; Ebel ing et aL]|1998| [Bohringer et al.||2001l |Mullis et al" 
2003 ) produced catalogs with hundreds of galaxy clus- 



ters, while pointed X- ray observations h ave discovered 
systems up to z ~ 1.4 (iMullis et aT|2005| ). 

With this article we inaugurate the Southern Cosmol- 
ogy Survey (SCS). This project, funded by the National 
Science Foundation under the Partnership in Research 
and Education (PIRE) program, is a multiwavelength 
(radio, mm-band, optical, UV, and X-ray) large area 
survey specifically coordinated with ACT observations 
of the southern sky. The goal of the SCS is to maximize 
the scientific return from the new ground-based CMB ex- 
periments and therefore focuses on specific observational 
studies relevant to this science, such as the selection func- 
tion of galaxy clusters across wavebands, cluster mass 
determination, and the establishment a "gold" sample 
of clusters for cosmology and galaxy evolution studies. 
Here we present results from an ~ 8 deg 2 optical imag- 
ing survey of the southern sky that overlaps the common 
SZ survey region. The purpose of this paper is twofold: 
(1) to present the details of our data reduction pipeline 
and analysis software and (2) identify new galaxy clus- 
ters, constrain their redshifts and masses, and predict 
their SZ signals. Photometric redshifts come from the 4- 
band imaging data, while our mass estimates are inferred 
from the optical luminosity (L200) and richness (iVfoo) °f 
the clusters, using relations calibrated by the Sloan Dig- 
ital Sky Survey (SDSS). For the eight massive clusters, 
out of 38 identified (37 are new sources) in the survey 
area, we present positions, richness estimates, masses, 
and predictions for the integrated Compton ^/-distortion 
of the SZE using empirical power-law relations based on 
iV-body simulations. Throughout this paper we assume 
a flat cosmology with Hq = lOOh km s _1 Mpc -1 , h = 0.7 
and matter density Q m = 0.3. 

2. DATASET AND METHODOLOGY 

Our study is based on the optical multi-band analy- 
sis of public data from the Blanco Cosmology Survey 9 
(BCS). This is a NOAO Large Survey Project that was 
awarded 45 nights over three years on the Blanco 4-m 
telescope at the Cerro Tololo InterAmerican Observatory 
(CTIO). The survey aims to image two 50 square-degree 
patches of the southern sky in four optical bands (griz) 
using the 8192 x 8192 pixel (0.36 deg 2 ) MOSAIC II cam- 
era in order to attain a sensitivity about an order of 
magnitude deeper than the SDSS imaging. The targeted 
areas are centered near declinations of —55° and —52° at 
right ascensions of 23 hr and 5 hr respectively; each of 
these patches lies within a larger common region of the 
southern sky that both ACT and SPT plan to survey. 
The BCS began in 2005 and has completed three years 
of data taking. For this paper we have processed and 
analyzed public data from the first year of the survey us- 
ing an independent software pipeline developed by us at 
Rutgers University. The data we present was obtained 
on 15 nights of observing near the end of November and 
the beginning of December 2005 and cover an area of 
c± 8 deg 2 in the 23 hr region. In Table [l] we show the ob- 
serving dates, photometric conditions, lunar illumination 
and observed bands for the 19 tiles that make up the full 
extent of the observations analyzed here. In the follow- 

9 http://cosmology.uiuc.edu/BCS/ 



TABLE 1 

2005 Observations in the 23hr Field 



# of Tiles Obs 



Date 


Photometric 


Lunar Ilium 


g 


r 


i 


z 


18 Nov 2005 


yes 


89.4% 


1.0 


1.0 


3.0 


3.0 


19 Nov 2005 


yes 


0.0% 


4.5 


4.5 


1.0 


1.0 


20 Nov 2005 


no 


0.0% 


8.0 


8.0 


0.0 


0.0 


22 Nov 2005 


no 


0.0% 


0.0 


0.0 


0.0 


0.0 


24 Nov 2005 


yes 


0.0% 


2.0 


2.0 


0.0 


0.0 


26 Nov 2005 


yes 


0.0% 


3.0 


3.0 


1.0 


1.0 


28 Nov 2005 


yes 


0.0% 


0.0 


0.0 


2.7 


2.7 


30 Nov 2005 


yes 


0.0% 


0.0 


0.0 


0.3 


0.3 


02 Dec 2005 


yes 


0.0% 


0.5 


0.5 


0.0 


0.0 


04 Dec 2005 


no 


15.3% 


0.0 


0.0 


0.0 


0.0 


05 Dec 2005 


no 


24.5% 


0.0 


0.0 


1.3 


1.7 


06 Dec 2005 


yes 


35.1% 


0.0 


0.0 


1.3 


1.0 


08 Dec 2005 


yes 


57.8% 


0.0 


0.0 


2.3 


2.7 


10 Dec 2005 


yes 


78.4% 


0.0 


0.0 


3.0 


3.0 


11 Dec 2005 


yes 


86.6% 


0.0 


0.0 


3.0 


3.3 



Note. — Observing conditions during the 2005 run of the BCS, 
consisting of only the 19 tiles that were fully observed in all 4 
bands in the 23hr region. Lunar illumination is the percentage at 
midnight local time in the direction toward the center of the region 
surveyed (R.A. 23 hours, decl. — 55.2deg). 



ing we describe the steps followed and tasks performed 
by the pipeline. 

2.1. The Rutgers Southern Cosmology Pipeline 

The Rutgers Southern Cosmology image analysis 
pipeline is written in Python with a scalable object- 
oriented design based on existing public astronomical 
software that is aimed at processing a large dataset in 
a repeatable, stable and semi-automated fashion. 

The initial standard image processing steps for e ach ob- 
servin g night are handled by the IRAF/mscred QValdes 
|1998[ ) procedures via the STScI/Pyraf interface. These 
include: overscan trim, bias correction, CCD cross-talk 
coefficients corrections as well as dome flat field cor- 
rection. The pipeline also executes secondary CCD 
calibration steps on the science images which include 
the creation of super sky- flats, fringe patterns for i— 
and z— bands and their corresponding correction and re- 
moval. Additionally, procedures affecting the cosmetic 
appearance of the images, such as cosmic ray rejection, 
removal of saturated star bleed- trails, and generation of 
bad pixel masks, are automatically performed at this 
stage. Astrometric re-calibration and WCS plate solu- 
tion are also handled automatically at this stage on the 
pre-stack science images using IRAF's mscred/mscmatch 
task by matching several hundred sources within each tile 
with stars from the US Naval Observatory Catalog. We 
achieve good accurate astrometric solutions (the residual 
error in matched source positions was typically < 0.1") 
as tested using the overlapping regions between neigh- 
boring pointings. Photometric standard star fields were 
processed together with the normal science images and 
photometric zero points for each observing night were ob- 
tained using a few hundred standards from the Southern 
Hemisphere Standards Stars Catalog ( |Smith et al.|2QQ7[ ). 
Like these authors we use the AB magnitude system. 

The survey strategy followed a predetermined observ- 
ing pattern, which typically consisted of exposures of 
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Fig. 1. — The mean lunar illumination fraction for each tile at 
time of observation in the 23hr field as a function of the i-band 
magnitude limit for 90% completeness. 



2 x 125s, 2 x 300s, 3 x 450s and 3 x 235s in the g, r, i 
and z— bands respectively with offsets of 3 — 5 arcmin 
(within each filter) intended to provide significant over- 
lap between neighboring MOSAIC II tiles and fill in the 
gaps between CCDs chips. We used the overlapping re- 
gions between tiles to adjust the photometric zero points 
of non-photometric nights using matched sources from 
adjacent photometric tiles. This ensured a homogeneous 
photometric calibration across the full survey region with 
typical variations below 0.02 mags. 

Image alignment, stacking and combination as well as 
catalog generation are performed at a secondary stage 
by the pipeline using association files, which describes a 
logical group of exposures and filters, created for each 
tile. Science images were mosaiced, al igned and median 



combined using SWarp 
of 0.266^/pixel. Source 



2006) to a plate scale 



Bert in 

etectioii and photometry mea- 
surements f or the science catalogs w ere performed using 
SExtractor ( |Bertin fc Arnouts|1996[ ) in dual-image mode 
in which sources were identified on the i— band images 
using a 1.5<r detection threshold, while magnitudes were 
extracted at matching locations from all 4 bands. 

As our data set is composed of observations taken over 
several weeks under varying conditions during the 2005 
campaign, we determined, for each tile individually, the 
i-band magnitude limit at which the galaxy detection 
limit was complete to 90%. To compute this limit we use 
the fact that the galaxy number counts follow a power- 
law function, which we fitted in the magnitude range 
19.5 < i < 21.5 in each tile and extrapolated to ob- 
tain the magnitude at which the galaxy number counts 
dropped by 10%. We took this as the 90% complete- 
ness limit for the tile. We found variations of roughly 1 
magnitude on the limits among the 19 tiles and, in an 
attempt to understand this, we investigated a possible 
correlation with lunar illumination at the time of obser- 



vation. Figure [T] shows a clear trend between the i-band 
magnitude limit and the lunar illumination. We report 
a mean limit i = 22.62 ± 0.25 and we set a conservative 
magnitude limit of i = 22.5 for our full catalog. 

2.2. Photometric Redshifts 
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Fig. 2. — Comparison between photometric and spectroscopic 
redshifts for 5000 galaxies in the SDSS DR6. The upper panel 
shows the comparison between spectroscopic and photometric red- 
shifts obtained under the maximum likelihood assumption. The 
lower panel uses a Bayesian assumption with a custom empirical 
prior on galaxy brightness for the photometric redshifts. Symbols 
are color coded according to the best determined SED by BPZ. 
SEDs for E/S0 galaxies tend to dominate at higher redshifts as 
the SDSS is biased towards early-type galaxies which are the most 
luminous population at these redshifts. 

^From the multi-band photometry the pipeline com- 
putes photometric redshifts and redshift probability dis- 
tributions Pbpz(z) for each object using the g,r,i,z 
isophotal magnitud es, as defined b y the i-band detection, 
and the BPZ code (penitez 20001). We use a magnitude- 
based empirical prior (Benitez, private communication) 
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taken from the SDSS and HDF-N spectroscopic red- 
shift distributions, which accounts for the tendency of 
fainter galaxies to be m ore likely foun d at higher red- 
shifts (See Fig 4. from |Bemtez| [2000 , as an example). 
Because the area covered by the available NO AO imaging 
does not include any publicly available spectroscopic red- 
shift information for z > 0.1 (NASA/IPAC Extragalac- 
tic Database, NED), we investigated the accuracy of our 
photometric redshift estimates using ancillary data. As 
the ability of BPZ to estimate photometric redshifts at 
fainter magnitudes (i > 20) from multi-band photometry 
has been consistently est ablished in the past using filter 
sets similar to ours here (Benitez et al. 1120041 iMobasher 



et al. 2004| Cross et al.||2004} for example), we focused 
on the redshift accuracy at z < 0.5. To this end we ex- 
tracted g, r, z, z photometry from the DR6 SDSS for 5000 
randomly-selected bright galaxies (r < 20 mag) with re- 
liable spectroscopic redshifts, matching the depth and 
signal-to-noise ratio of our galaxy sample. We computed 
photometric redshifts for the SDSS spectroscopic sample 
in the same way as just described and compared the re- 
sulting values to the spectroscopic redshifts. We found, 
not surprisingly, that simply employing the maximum- 
likelihood (ML) condition is an ill-suited approach for 
redshifts below z < 0.3 in the absence of a bandpass bluer 
than 3000A as it largely over-estimates redshifts and pro- 
duces an unacceptable number of catastrophic outliers. 
Recently, Niemack et al. (2008) have demonstrated how 
the addition of bluer bands using GALEX UV imaging 
greatly improves ML estimates and reduces the need for 
priors. On the other hand, Bayesian estimates give re- 
sults with typical rms errors of 5z ~ 0.02 and with almost 
no catastrophic outliers. In Fig. [2] we show the results 
of our comparison between ML and Bayesian photomet- 
ric versus spectroscopic redshifts color-coded according 
to the spectral energy distribution (SED) determined 
by the BPZ code. It is clear from the figure that at 
higher redshifts the SDSS population is dominated by 
early types as these tend to be the most luminous ob- 
jects. We also note (see fig. pi lower panel) that on av- 
erage our photometric redshnts tend to overpredict the 
true redshifts. The mean bias level (z spec — ^bpz) for 
galaxies with E/SOs SEDs is largest around z^ 

ec 00 0-3 

where it is on the order of Sz ~ —0.03. In table|2]we show 
the mean bias and standard deviation for three redshift 
intervals for z spe c — z B pz as well as the standard dz de- 
fined as dz = z spec — ^BPz/(l + ^spec)- In summary we are 
able to determine the redshifts for early-type galaxies to 
an accuracy better than 0.1 across the redshift range of 
the survey. This is encouraging since early-type galaxies 
are the predominant population in clusters of galaxies 
and good photometric redshift determination is essential 
for successful cluster finding, as we discuss in the next 
section. 

In Figure [3] we show the photometric redshift distribu- 
tion for all galaxies within our flux completeness limit, 
i < 22.5, as well as the filter responses for the survey 10 . 
Our distribution peaks around z ~ 0.6, which sets a con- 
servative upper limit to the redshift at which we are able 
to detect optical clusters. 

2.3. Computing Overdensities and Finding Clusters 

10 http:/ /www. ctio.noao.edu/instruments/FILTERS/index. html 
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Fig. 3. — The photometric redshift distribution of the galaxy 
sample used for finding clusters (i.e.: galaxies i < 22.5 mag). 



TABLE 2 

Photometric Redshift Simulations Statistics 



Redshift 


(Zspec — ZBPZ) 


(Zspec — ZBPz)rms 


(dz) 


cr z 


0.0 - 0.2 


-0.017 


0.042 


-0.015 


0.038 


0.2 - 0.4 


-0.027 


0.059 


-0.020 


0.047 


0.4-0.6 


-0.002 


0.070 


-0.001 


0.048 



Note. — The mean difference and standard deviation between 
spectroscopic and recovered photometric redshifts as well as for dz 
for all galaxies with SED determined to be E/SOs in three redshift 
ranges. 



One of the main goals of the current SZE experiments 
is to define a mass-selected sample of galaxy clusters out 
to large redshifts. At long l ast this is beginning to happen 
( Staniszewski et al. 2008; Me nanteau fc Hughes| |2009), 
after a number of successful individual detection s of the 
SZE in well-known optical or x-ray clusters (see |Birkin 



shaw et al. 



Reese et al. 



1991; Jo nes et al.|p93j |Herbig et ah} 1995; 
2000 ; |G6mez et a l. 2003; Bonam ente et al. 



2006| Halverson et al. |2008| and references therein). If 
we want to understand the systematics of SZE surveys 
it is crucial to compare with cluster identifications using 
independent methods. In this section, we describe our ef- 
fort to select clusters of galaxies from multi-wavelength 
optical imaging. There are several methodologies and 



a plethora of papers describing these te c hniques ( Post- 
man et al.||1996l |G ladders fc Yee||2005l |Koester et aT 
2007} |Eisenhardt et al ||2008[ for example) but they all 
rely on the same well-known properties of galaxy clus- 
ters: a) early-type galaxies are the dominant population, 
b) cluster galaxies have very similar colors, and display 
tight color-magnitude relationships across several orders 
of magnitude in luminosity, and c) the surface number 
density of cluster galaxies falls off with distance from the 
center roughly as a power law P(r) oc l/r a . We search 
for clusters using a matched filter appro ach similar to the 
one described in |Postman et al. ( 1996[ ) and then define 
membership and estimate richness of the clust ers using 
the MaxBCG prescription (Ko ester et al.|[2007| ). 

Our cluster finder method folds in the contributions 
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from a cluster spatial profile filter function P(r), a lu- 
minosity weight L(m) and the BPZ redshift probability 
distribution Pbpz(^) from each source to generate likeli- 
hood density maps (at pixel positions denoted by z, j) or 
a "filtered" galaxy catalog S(i,j)(z) over the area cov- 
ered by the survey as a function of redshift, namely, 



S(iJ)(z) = Y^P(r k [iJ]) L(m k ) f + ^ p B pz( 

k=l Jz-Az 



(z k )dz. 



Specifically we use a profile with the form 

/ / \ 1 1 

P(r/r c ) 



- 0, otherwise 

which is normalized as 

P(r/r c )27rrdr = 1 



/ 

Jo 



(i) 

, if r < r cut 
(2) 
(3) 



and where r c is the typical cluster core radius and r cut is 
the cutoff limit for the function. In our analysis we chose 
r c = 175 kpc and r cut = 10r c . We also use a luminosity 
weight L(m) given by 



L(m) = 



<j)(m 



)10 -0.4(m-m*) $( m . 



b(m) 



b(m) 



(4) 



where m* is the apparent magnitude corresponding to 
M*. This function is normalized as 

<£(m — m*)dm — 1 (5) 



f 

Jo 



where mu m is the flux limit of the sample (i = 22.5), 
b( m) is the number of background galaxies and <p(m) is 
the |Schechter| ( |1976[ ) galaxy lum inosity function. We use 
the parameters computed from |Brown et al.| ( |2QQ7[ ) for 
the evolving luminosity function of red galaxies, with a 
faint-end slope a = —0.5 and M*(z) between < z < 1. 
For our estimation of bim) we use the number counts 
from Yasuda et al. (2001 ). We generate likelihood density 



maps with a constant pixel scale of 1.2 arcmin at Az = 
0.1 intervals between 0.1 < z < 0.8 over the surveyed 
regions. In Figure [4] we show an example of a likelihood 
density map centered at z = 0.2 on which we superpose 
outlines of the 19 tiles that define the region studied here. 

Cluster candidates are selected from the peaks of the 
likelihood density maps. In order to define peaks consis- 
tently we constructed noise maps by randomizing the po- 
sitions of the input catalog to produce likelihood density 
maps following exactly the same procedure as described 
above. These maps represent the noise floor level above 
which we desire to detect clusters. We define our initial 
list of cluster candidates from 2a peaks in the likelihood 
density maps, where a is defined as the median value 
in the noise maps. Cluster candidates were checked for 
multiple detection in different likelihood maps. To avoid 
duplication we considered a system unique if detected 
in two adjacent redshift maps and with the same center 
(i.e., within 3'). 

2.4. Contamination and Completeness 
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Fig. 4. — The computed likelihood density map image centered 
at z = 0.2 and width Az = 0.1 over the 23hr field. Bright regions 
in the image represent denser areas. The red lines represent the 
area covered by each of the 19 tiles that comprise the area studied. 



We perfomed simulations to investigate our cluster se- 
lection function by estimating the contamination and re- 
covery rates of our cluster finding technique. The lack of 
distance information in imaging surveys is the principal 
source of contamination as fluctuations in the projected 
two-dimensional galaxy distribution as well as random 
alignments of poor groups may result in false apparent 
overdensities. We explor e this issue following the s ame 
methodology as used by |Postman et ah] ([2002 ) and |Gal| 
et al.| ( |2003[ ) which rely on generating Monte Carlo rep- 
resentations of the galaxy sample with an angular t wo- 
point corr elation function si milar to that observed (In- 
fante||1994| ). As described in |Postman et al.| ( |2002[ ) we 
implemented the Ray leigh- Levy (RL) random walk pro- 
cess using Mandelbrot's (1975) elegant fractal prescrip- 
tion to simulate galaxy positions on the sky, such that 
galaxy pairs are placed in a randomly chosen direction 
at distance 6 drawn from the distribution: 



if > O 
if < O) 



(6) 



where we chose 6q and d to match the observed galaxy 
distribution of our sample. In practice we generate simu- 
lated distributions by starting from a randomly selected 
location within the survey boundaries and generate posi- 
tions following the RL random walk allowing up to seven 
galaxies to be drawn around this location. We then se- 
lect a new center randomly and the process is repeated 
until we generate the same number of galaxies as in the 
observed sample. We then process the RL distribution to 
generate likelihood maps using the same procedure and 
parameters as for the real data and use these to investi- 
gate the rate of false detections as a function of estimated 
redshift. Since the RL distributions by construction do 
not explicitly include clusters, we assess the false posi- 
tive cluster detection fraction by taking the ratio of de- 
tections in the simulations per area unit to the observed 
number of candidates in the real data. We find that at 
low redshift the false positive fraction is zero (there is 
virtually no contamination), while at redshifts of z = 0.6 
and z = 0.7 the false positive fraction grows to values of 
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Fig. 5. — The cluster recovery fraction as a function of red- 
shift as extracted from the simulations for clusters with N s i m = 
15, 20, 25, 35, 50, 80 and 120 galaxies and a profile with slope a = 
1.8 and core radius 0.150 Mpc. 



1% and 19%, respectively. We conclude that false posi- 
tives are not an important source of spurious detections. 

We investigate the selection function for our galaxy 
sample by simulating galaxy clusters of various richness 
and shapes at different redshifts and examining their re- 
covery fractions. Specifically we generated clusters with 
random ellipticities uniformly between 0.1 < z < 0.7 
using an r~ a pro file for the galaxy distribution (|Lubin| 



Postman| |1996) with a = 1.8 and r c = 0.150 Mpc 
and using th e lumi nosity function for red galaxies from 



Brown et al. (2007). These clusters are inserted 20 at a 
time in the observed catalogs at random positions and 
redshifts, but avoiding the locations where clusters were 
detected. In total we generate 10,000 simulated clusters 
with richness values, N sirn , of 15, 20, 25, 35, 50, 80 and 
120 galaxies uniformly distributed in redshift. 11 We pro- 
cess each realization using the same parameters as for 
the observed data and record the number of clusters re- 
covered as a function of redshift and galaxy richness. In 
Figure [5] we show the results of this exercise where we 
plot the recovery fraction as a function of redshift for 
the seven cluster richnesses simulated. We conclude that 
for the rich clusters N sim > 50 we are always nearly com- 
plete (80 — 90%) for z < 0.6 while for poorer clusters we 
only detect at best ~ 30% around z = 0.3. 

11 These richness values fold in the flux limit of the survey and 
the membership prescription as described in the next section, so 
that they are roughly comparable to the N ga i values we give for 
the detected clusters. 



3. CLUSTER PROPERTIES 

One of our main drivers in searching for clusters at op- 
tical wavelengths is to correlate them with SZ detections 
in the new blind SZ surveys. The signal to be detected in 
the mm-band experiments (i.e., the ^/-distortion due to 
inverse Compton scattering) is related to the number of 
hot electrons in t he intracluster medium ,and simulations 
have shown (see|Motl eta l. 2005; Nagai||2006| ISehgal et| 
al.|2007 Bhattacharya et al. 2 007[ for example) that the 
SZE signal is closely linked to cluster mass. Our analysis 
of the optical survey has provided positions, redshifts, 
and fluxes of galaxies, from which we infer the under- 
lying cluster mass using scaling relations established by 
the SDSS, where cluster masses were determined from 
weak lensing. 

3.1. Defining Cluster Membership 
The current state-of-the-art mass tracers for clusters 



of galaxies using optically observe d parameters (John- 



ston et al.|[2007 



Reyes et al. 2008) have been extracted 



from a sample of around 13, 000 opticall y-selected clus- 
ters from the SDSS MaxBCG catalog (iKoester et al. 



2007). In this paper, we apply these scaling relations 



to our cluster sample and obtain mass estimates from 
which we additionally predict SZ distortions. To be fully 
consistent with previous analyses, we define membership 
and all other rel evant cluster observ able s following the 



Koester et al. 



same method as |Reyes et al.| ( |2008| ) and 
( 2007| ). 

We begin by examining each candidate-cluster peak 
in the density maps and select the brightest elliptical 
galaxy in the cluster (BCG), which is taken to be the 
initial center and redshift z Q of the system. We then 
use galaxies defined as E or E/SOs (i.e., BPZ SED types 
and 1 only) within a projected radius of 0.5/i _1 Mpc 
and redshift interval \z — z Q \ = \Az\ = 0.05 to obtain 
a local color-magnitude relation (CMR) for each color 
combination, g — r, r — i and i — z, and the cluster 
mean redshift, z c , for all cluster members, using a 3<r 
median sigma-clipping algorithm. We use these to deter- 
mine A^Mpc? the number of galaxies within l/i _1 Mpc of 
the cluster center. Formally, we compute N ga i = A^im pc 
by including those galaxies within a projected lh~ 1 Mpc 
from the cluster center that satisfy three conditions: (a) 
the galaxy must have the SED of an E, E/S0 according 
to BPZ, (b) it must have the appropriate color to be a 
cluster member (i.e., colors within 3a of the local CMR 
for all color combinations) and (c) it must have the right 
luminosity (i.e., dimmer than the BCG and brighter than 
0.4L*, where we use the correspon ding absolute magni- 



(20C 

bp; 



2007) redshifted with the 



FZj . We designated clus- 



tude M* from Brown et al. 
elliptical SED template from 

ter members according to the estimated cluster size i?200> 
defined as the radius at which the cluster galaxy density 
is 200^~ 1 times the mean space density of galaxies in the 
present Universe. We estimated t he scaled radius R20 



= O.lbQN^^h^Mpc whic 



^200 = 0.156N\*i®^h x Mpc which is derived from the 



using the empirical relation from [Hansen et al.| (|2005j) 

ich 



SDSS and we assume it holds beyond z ~ 0.3 for our 
higher redshift clusters. 

In our analysis we use AT|qq, L200, and Lbcg to scale 
clus ter op tical parameters with mass, following Reyes et| 



al. (120081). The cluster richness, iVfoo, is the number 
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Fig. 6. — The composed gri color image (left panel) and color-magnitude relations (right) for X-ray cluster RXJ2325.6— 5443 from[Mullis 
|et al.| ( |2b03|). Red points represent galaxies classified as E/S0 by BPZ that satisfy the conditions to be cluster members as described in the 
main text. Black dots are non-member galaxies in a 5 7 region near the cluster center. 



of E/S0 galaxies within i? 2 oo with colors and luminosi- 
ties that satisfy conditions (b) and (c) above. Similarly, 
L200 is the total rest-frame integrated r-band luminos- 
ity of all member galaxies included in ATfoo m units of 
lO lo /z~ 2 L and Lbcg is the rest-frame r-band luminos- 
ity of the BCG. 

In order to have reliable estimates of 7V gal it is neces- 
sary to determine the galaxy background contamination 
and implement an appropriate background subtraction 
method. The lack of spectroscopic redshifts in our sam- 
ple only allowed a statistical removal of unrelated field 
galaxies with similar colors and redshifts that were pro- 
jected along the line of sight to each cluster. We as- 
sumed that the presence of a cluster at some redshift is 
independent of the field population seen in projection. 
Therefore we estimate the surface number density of el- 
lipticals in an annulus surrounding the cluster (within 
^200 < r < 2i?20o) with Az = 0.05 and the same colors 
as the cluster members. We measure this background 
contribution around the outskirts of each cluster and ob- 
tain a corrected value N ga i which is used to compute R200 

and then corresponding values of iVfoo an d £200- The 
magnitude of the correction ranges between 15—20%. We 
will refer to the corrected values hereafter. Moreover as 
our analysis is based on a magnitude-limited sample it is 
worth considering the fraction of lower luminosity galax- 
ies that will fall below our magnitude limit (i = 22.5) at 
higher redshifts. As in Menanteau & Hughes (2009), if 



we make the assumption that the cluster population is 
like that of the five clusters at z < 0.2 in our sample (see 
Table [5| and M* evolves passively, then we can compute 
the fraction of L 2 oo missed for clusters at higher red- 
shifts. We estimate that we are missing 4%, 13%, 29% 
and 38% of the cluster luminosity at z = 0.4, 0.5, 0.6 and 
0.7 respectively. Given the uncertainty in this correction 
factor, we do not include it in our quoted luminosity 
values for the higher redshift clusters. This means our 
cluster masses, M(L20o), are underestimated by roughly 
these factors. 

3.2. Recovery of Known Clusters 

The area covered by the BCS in the 23hr region is a 
virtual desert in terms of known clusters and spectro- 
scopic redshifts for galaxies with z > 0.1. We found 
one catalogued X-ray selected clus ter from the 160 deg 2 
ROSAT survey flMullis et al. 120031): RXJ2325.6-5443 at 



0.102 with an X-ray flux of F x = 2.2 x 10~ 13 erg 
cm -2 s _1 in the 0.5-2 keV band. Our cluster finding al- 
gorithm easily recovered this cluster and produced a pho- 
tometric redshift estimate of (z) = 0.10 ± 0.02. Figure [6] 
shows the gri color composite optical image of the cluster 
as well as the color magnitude diagrams for cluster mem- 
bers. Using the techniques described below we estimate 
the mass of RXJ2325.6-5443 to be M(L 200 ) = 2.1 x 10 14 
which is just below the detect ability limit of ACT and 
therefore will not be included in our SZE predictions. 
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Fig. 7. — The mass observable parameters for the 38 clusters in our sample. Filled circles represent the eight new clusters above the 
mass limit M(L2oo) > 3 x 10 14 Mq. The left panel shows the relation between the mass predictions using L200 vs. N^qq as the main 
parameter; the agreement is generally to within a factor of two, except for the nearest systems. The central and right panel show the mass 
estimates as a function of the cluster redshift using L200 and N^qq respectively, in units of solar mass. 



From the M-T x (|Evrard et al.||1996[ ) and L x -T x ( [Ar^ 
naud fc Evrard||1999p relations we estimate a tempera- 
ture oi kT ~ 2.2 keV and an X-ray flux of F x ~ 5 x 10" 13 
erg cm -2 s _1 (0.5-2 keV band) which is in rough agree- 
ment (factor of 2) with the published value. Considering 
the entire sample of 38 clusters, RXJ2325.6— 5443 is the 
closest and, based on our estimated masses from the op- 
tical properties and the M-T x and L x -T x relations, has 
the highest predicted X-ray flux of the sample. Still, even 
this value is below the X-ray detection threshold of the 
ROSAT All Sky Survey (RASS) and so, as expected, we 
find that none of our new optical clusters are significant 
RASS X-ray sources. 

3.3. Cluster Mass Estimation 

In this section we use the mass-richness relations based 
on A/impc to wei gh o ur new optica l clust ers. Both|John^ 



ston et al.| fl2007D and |Reyes etaL] ( |2008| ) found that the 
luminosity-mass and richness-mass relations were well 
described by power-law functions and they measured the 
normalizations and slopes in these relations using y 2 
minimization. Their values are in broad agreemen t, but 
we will use the fits provided by Reyes et al. ( 2008| ) since 
they restrict their fits to clusters with A^oo > 10 and 
give results for two redshift bins: 0.10 < z < 0.23 and 
0.23 < z < 0.30. It is important to note that our clus- 
ters go well beyond z — 0.3 and that we extrapolate the 
relation in the last redshift bin for clusters with z > 0.3. 
We investigated the two fitting fu nctions based on L 200 

for 



2008 



and A^Iqq, (see section 5.2.1 from |Reyes et al. 
full details), which are described as: 



M(7V 200 ,L B cg) = M^(iV 200 /20) a -(L B cG/4cG) 7Ar 



M(L 200 ,£bcg) = M£(L 20 o/40)^(L B cg/4cg 



\1L 



(7) 
(8) 



where M is the mass observational equivalent of M 2 oop 12 
in units of 1O 14 M , L200 is in units of lO lo /i _2 L and 
the I/bcg dependence is normalized by its mean value. 

12 M200P is the halo mass enclosed within R200, defined as a 
radius of spherical volume within which the mean density is 200 
times the critical density. 
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Fig. 8. — The gray circles show the total number of clusters with 
M(L2oo) > 3 x 10 l4 Mq observed in the ~ 8 deg 2 region analyzed 
here. The thin and thick lines show the predicted number of clus- 
ters in the region for cluster masses greater than 2 x 10 14 Mq and 
3 x 10 14 Mq respecti vely, using the WMA P 5-year best-fit cosmo- 
logical model and the [Jenkins et al.| ( |2001| ) halo mass function. 



This is also described by a power-law function for a given 



value of L200 and iVjgJ: 



J BCG — 



^BCG = ^BCG(^20o) = a L^200 



'bcg(JV2oo) = ajvJVaoo ( 9 ) 

(10) 

The published best-fitting parameters for M°, a and 7 in 
Eqs. ([7| and (8| as well as the values of a, b for Eqs. (|9| 
and ([l0 ) are shown in Table [3] The combination of equa- 
tions^), <Js]) and |9]), ([Io| for JVjgJ and L200 respectively 
enable us to obtain mass estimates for any cluster with 
AUi > 10. 
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TABLE 3 

Mass-richness power-law function Best fitting parameters 



Redshift 


(io 10 /i- 

a N 


- 2 L Q ) 

a L 


b N 


b L 


K 


a N 


IN 


M° 


a L 


IL 


0.10 < z < 0.23 
0.23 < z < 0.70 


1.54 
1.64 


7.77 
7.92 


0.41 
0.43 


0.67 
0.66 


1.27 ±0.08 
1.57 ±0.14 


1.20 ±0.09 
1.12 ±0.15 


0.71 ±0.14 
0.34 ±0.24 


1.81 ±0.15 
1.76 ±0.22 


1.27 ±0.17 
1.30 ±0.29 


0.40 ±0.23 
0.26 ±0.41 



The left panel of Fig. [7] compares the masses obtained 
using A^2oo and £200 f° r °ur 38 optical clusters. The solid 
line denotes equality between the estimates while the two 
dashed lines show a factor of 2 range. Most clusters fall 
within this range, establishing a lower bound on our mass 
error. The two most significant outliers correspond to 
nearby clusters (see the right panel of Fig. [7| for which 
7V200 is evidently overpredicting the mass compared to 
£200- I n one °f these cases (which is RXJ2325.6— 5443), 
the mass inferred by A^2oo grossly overpredicts (by more 
than an order of magnitude) the estimated X-ray flux. 
The natural concentration of luminosities in clusters (i.e., 
that more luminous galaxies dominate the central regions 
of clusters) suggests that the mass estimate derived from 
L200 should be more robust than that from A^oo as a 
function of redshift. For these reasons we use the cluster 
mass estimate derived from L 2 oo m predicting the SZE 
signal. Hereafter we refer to this mass as M^qq. 

Table [4] displays the 8 clusters with M(L 2 oo) > 3 x 
1O 14 M in the 8 square-degree sky area covered by our 
analysis. These clusters also have M(7V 2 oo) above the 
same mass limit. In Figure [8] we show the cumulative 
number distribution of these clusters as a function of 
redshift, compared with expectations from simulations 
(H. Trac, private communication) , using the mas s func- 
tion of dark matter haloe s from i Jenkins et aL] ( |20Q1[ ) 
and WMAP5 cosmology ( |Komatsu~et al.||2008f . Tne 
heavy solid curve shows expectations for a mass limit 
of 3 x 1O 14 M , while the light solid curve indicates a 
mass limit of 2 x 1O 14 M . This figure demonstrates that 
the number of clusters we observe is consistent with a 
mass limit in the range of 2 - 3 x 10 14 M o . In Table [5] 
we display the properties for the remaining clusters witn 
M(L 20 o) < 3 x 1O 14 M . 

4. PREDICTIONS OF THE SUNYAEV-ZEL'DOVICH 
EFFECT SIGNAL 

Scaling relations between the integrated thermal SZE 
signal and cluster mass have emerged from current N- 
body plus hydrodynamic simulations of galaxy clusters. 
In this section we use these relations and our optically- 
derived mass estimates, MJqq, to predict the SZE signal 
to be observed by ACT and SPT when these experiments 
survey this sky region. 

The SZE signal consists of small distortions to the 
CMB spectrum originating from inverse Compton scat- 
ter ing_by_d^c_tron^jii_^^ of clusters of galax- 
ies QSunyaev fc Zeldovich||1980| ). Here we consider the 
thermal SZ flux K, defined as the integrated Compton 
parameter, 

Y = d\{z) f y<m= k ^r f n e T e dV, (11) 

JQ Jv 

where n e and T e are the number density and temperature 



of hot electrons in the cluster, m e is electron rest mass, c 
is the speed of light, <jt is the Thompson scattering cross 
section and (Ia(z) is the angular diameter distance. The 
projected area, cL4, and solid angle, dQ, are related by 
the angular diameter distance as dA = d 2 A (z )d£t. 

Self-similar scaling relations ( |Kaiser|1986[ ) predict that 
the virialized mass M in clusters scales with the gas tem- 
perature as M oc T 3 / 2 /E(z) where E(z) = (O m (l + 
z) 3 + ^a) 1 ^ 2 for a flat cosmology. If clusters were 
isothermal, we would expect their SZE signal to scale 
like Y oc /gas^haio^ and therefore the self-similar SZ 
flux-mass scaling relation should have the shape Y oc 
f gas M 5 / 3 E 2 / 3 (z) where / gas is the cluster mass fraction. 
However, clusters are not always isothermal or in hydro- 
static equilibrium and physical processes like star forma- 
tion and feedback will also contribute to deviations from 
self-similarity. 

Prompted by the upcoming SZ surveys, several studies 
have characterized in detail the Y — M scaling relation 
using simple power-law fits to cosmological N-body sim- 
ulations. Here we will use the recent Y — M fits from 



Sehgal et al.| (2007), who included gas simulations em- 
ploying small-scale cluster physics such as star formation 
and feedback into a large cosmological N-body simula- 
tion. ^From their catalog o f ~ 10 5 simu l ated clusters 
with M 200 > 7.5 x 1O 13 M , |Sehgal et ah] ( |2007| ) fit the 
relation 

^200 f ^200 



E(z) 2 / 3 



Vl0 14 M G 



(12) 



where Y200 and M200 are the projected SZ Compton 
^/-parameter and mass respectively in a disk of radius 
^200- We use the best fit values a = 1.876 ± 0.005 
and (3 = -5.4774 ± 0.0009 for all clusters regardless of 
redshift, as the redshif t dependence of th e fits is very 
weak (see Table 2 from Sehgal et al. 2007). The power- 
law index that they report is slightly steeper than ones 
quoted by some previous hydrodynamic simulations (see 



Hernand e~Monteagudo et al.|2006j [N agai 2006; M qtl et 
|al.l |2005[ for example), which |Sehgal et al.| ( |2007| ) at- 
tribute to their more realistic feedback prescription (i.e., 
including the effects of AGN and supernovae). We use 
the power-law model from equation [12] and the optical 
Mjoo mass estimates for our 8 new massive clusters to 
predict the integrated Y200 signal to be observed by ACT 
and SPT. Table [4] gives the results in physical units of 
Mpc 2 and observable ones of arcmin 2 . 

5. RESULTS AND CONCLUSIONS 

In this article we have laid out the techniques and 
methods for our analysis of a large multi-band optical 
survey with the Blanco telescope and Mosaic-II instru- 
ment under the aegis of the Southern Cosmology Survey. 
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TABLE 4 

Optical Clusters with M(L 2 oo) > 3 x 1O 14 M 











iV 200 




[M Q ] 


[Mpc 2 ] 




[arcmin 


2 ] 


ID 




z 




£200 [Lq] 


M(N 20 o) 


M(L 20 o) 


Y 20 o/E(z) 


2/3 


^200 




SCSO J233430.2- 


-543647.5 


0.35 


32.1 ±5.8 


43.8 ±7.0 


1.8 x 10 12 ± 2.9 x 10 10 


6.1 x 10 14 


3.6 x 10 14 


3.5 x 10" 


-05 


4.8 x 10" 


-04 


SCSO J233556.8- 


-560602.3 


0.52 


31.4 ±5.9 


33.5 ±6.4 


3.3 x 10 12 ± 1.4 x 10 11 


4.6 x 10 14 


7.2 x 10 14 


1.2 x 10" 


-04 


1.1 x 10" 


-03 


SCSO J233425.6- 


-542718.0 


0.53 


26.1 ±5.4 


26.9 ±5.6 


1.8 x 10 12 ±9.5 x 10 10 


3.4 x 10 14 


3.4 x 10 14 


3.2 x 10" 


-05 


3.0 x 10" 


-04 


SCSO J232211.0- 


-561847.4 


0.61 


31.1 ±5.7 


34.1 ±6.1 


2.8 x 10 12 ± 1.1 x 10 11 


4.6 x 10 14 


5.6 x 10 14 


7.9 x 10" 


-05 


6.8 x 10" 


-04 


SCSO J233731.7- 


-560427.9 


0.61 


25.6 ±5.4 


23.1 ±5.3 


1.7 x 10 12 ± 1.7 x 10 11 


3.0 x 10 14 


3.0 x 10 14 


2.6 x 10" 


-05 


2.2 x 10" 


-04 


SCSO J234012.6- 


-541907.2 


0.62 


22.7 ± 5.2 


24.9 ±5.4 


2.4 x 10 12 ± 1.4 x 10 11 


3.8 x 10 14 


5.2 x 10 14 


6.9 x 10" 


-05 


5.8 x 10" 


-04 


SCSO J234004.9- 


-544444.8 


0.66 


30.7 ±5.8 


37.5 ±6.7 


2.3 x 10 12 ±2.4 x 10 11 


4.3 x 10 14 


4.2 x 10 14 


4.6 x 10" 


-05 


3.8 x 10" 


-04 


SCSO J232829.7- 


-544255.4 


0.68 


74.6 ±8.7 


65.7 ±9.7 


4.4 x 10 12 ±4.3 x 10 11 


7.1 x 10 14 


8.3 x 10 14 


1.6 x 10" 


-04 


1.3 x 10" 


-03 



Note. — Catalog of the optical clusters with mass estimates > 3 x 1O 14 M0 from the M(L2oo) values. Each cluster's redshift is the mean 
photometric redshift computed using the elliptical in the center of the cluster. The ID is based on the position of the BCG. 



TABLE 5 

Optical Clusters with M(L 2 oo) < 3 x 1O 14 M 



[M ] 





ID 










z 






iV 200 










£200 [Lq] 






M(N 20 o) 


M(L 


200) 


SCSO 


J232540 


2 


-544430 


9 





10 


123.9 ±11 


2 


278.1 ± 18 


5 


3 


3 


X 


10 12 


±2 


8 


X 


10 10 


8 


6 


X 


10 14 


2 


1 


X 


10 14 


SCSO 


J232230 


9 


-541608 


3 





12 


69.5 ±8 


4 


145.3 ± 12 


8 


1 


7 


X 


10 12 


± 1 


1 


X 


10 10 


1 





X 


10 15 


1 


6 


X 


10 14 


SCSO 


J233000 


4 


-543707 


7 





14 


38.0 ±6 


4 


39.0 ±7 


2 


1 





X 


10 12 


± 1 





X 


10 10 


4 


3 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J232419 


6 


-552548 


9 





18 


39.3 ±6 


5 


36.9 ±7 





1 


4 


X 


10 12 


± 1 


8 


X 


10 10 


2 


5 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J233106 


9 


-555119 


5 





19 


26.3 ±5 


5 


35.3 ±6 


4 


9 


7 


X 


10 11 


±2 


9 


X 


10 10 


1 


1 


X 


10 14 


5 


5 


X 


10 13 


SCSO 


J233252 


9 


-561454 


1 





20 


34.0 ± 6 


1 


43.4 ±7 


2 


1 


5 


X 


10 12 


±2 


2 


X 


10 10 


2 


2 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J233215 


5 


-544211 


6 





20 


43.5 ±6 


8 


42.6 ±7 


6 


1 


8 


X 


10 12 


±3 


4 


X 


10 10 


3 


3 


X 


10 14 


1 


7 


X 


10 14 


SCSO 


J233037 


1 


-554338 


8 





20 


27.5 ±5 


6 


35.2 ±6 


4 


1 


1 


X 


10 12 


± 1 


7 


X 


10 10 


2 


7 


X 


10 14 


9 


9 


X 


10 13 


SCSO 


J233550 


6 


-552820 


4 





22 


14.5 ±4 


4 


10.6 ±3 


5 


7 


4 


X 


10 11 


±2 


2 


X 


10 10 


8 


3 


X 


10 13 


6 


6 


X 


10 13 


SCSO 


J232200 


4 


-544459 


7 





27 


34.6 ± 6 





41.0 ±6 


9 


1 


2 


X 


10 12 


± 1 


8 


X 


10 10 


3 


9 


X 


10 14 


1 


7 


X 


10 14 


SCSO 


J233522 


6 


-553237 








29 


31.4 ±5 


9 


32.1 ±6 


3 


1 


5 


X 


10 12 


±2 


4 


X 


10 10 


3 


2 


X 


10 14 


2 


2 


X 


10 14 


SCSO 


J233807 


5- 


-560304 


9 





30 


32.0 ±5 


9 


37.7 ±6 


6 


1 


6 


X 


10 12 


±3 


2 


X 


10 10 


4 


2 


X 


10 14 


2 


6 


X 


10 14 


SCSO 


J232956 





-560808 


3 





32 


39.6 ±6 


5 


37.0 ±6 


7 


1 


3 


X 


10 12 


±2 


3 


X 


10 10 


3 


9 


X 


10 14 


2 





X 


10 14 


SCSO 


J232839 


5- 


-551353 


8 





32 


40.3 ±6 


5 


18.9 ±5 


2 


7 


9 


X 


10 11 


±2 


2 


X 


10 10 


1 


7 


X 


10 14 


1 





X 


10 14 


SCSO 


J232633 


6 


-550111 


5 





32 


74.3 ± 8 


7 


35.2 ±7 


4 


1 


9 


X 


10 12 


±4 


5 


X 


10 10 


3 


2 


X 


10 14 


2 


8 


X 


10 14 


SCSO 


J233753 


8- 


-561147 


6 





33 


33.2 ±5 


9 


41.3 ±6 


7 


1 


8 


X 


10 12 


±3 


3 


X 


10 10 


4 


2 


X 


10 14 


2 


9 


X 


10 14 


SCSO 


J232156 


4 


-541428 


8 





33 


20.1 ±4 


8 


19.9 ±4 


7 


8 


6 


X 


10 11 


± 1 





X 


10 10 


2 


1 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J233003 


6 


-541426 


7 





33 


29.6 ±5 


7 


30.4 ±5 


9 


6 


6 


X 


10 11 


± 1 


8 


X 


10 10 


2 


9 


X 


10 14 


8 


8 


X 


10 13 


SCSO 


J233231 


4 


-540135 


8 





33 


45.9 ±6 


9 


42.7 ±7 


1 


1 


2 


X 


10 12 


±2 


2 


X 


10 10 


4 


1 


X 


10 14 


1 


7 


X 


10 14 


SCSO 


J233110 


6 


-555213 


5 





39 


21.1 ±4 


9 


20.8 ±4 


9 


7 


3 


X 


10 11 


± 1 


7 


X 


10 10 


2 


2 


X 


10 14 


1 





X 


10 14 


SCSO 


J233618 


3 


-555440 


3 





49 


17.4 ±4 


6 


15.3 ±4 


2 


6 


3 


X 


10 11 


±3 


6 


X 


10 10 


1 


8 


X 


10 14 


9 


4 


X 


10 13 


SCSO 


J233706 


3 


-541903 


8 





51 


25.5 ±5 


3 


29.9 ±5 


9 


1 


2 


X 


10 12 


±8 


5 


X 


10 10 


2 


6 


X 


10 14 


1 


6 


X 


10 14 


SCSO 


J233816 


9 


-555331 


1 





52 


19.8 ±4 


7 


19.8 ±4 


7 


9 


7 


X 


10 11 


±3 


4 


X 


10 10 


1 


8 


X 


10 14 


1 


3 


X 


10 14 


SCSO 


J232619 


8- 


-552308 


8 





52 


18.8 ±4 


7 


18.2 ±4 


5 


8 


1 


X 


10 11 


±6 


1 


X 


10 10 


2 


1 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J232215 


9 


-555045 


6 





56 


11.0 ±4 





7.4 ±2 


9 


1 


2 


X 


10 12 


±6 


1 


X 


10 10 


1 


1 


X 


10 14 


2 


4 


X 


10 14 


SCSO 


J232247 


6 


-541110 


1 





57 


14.9 ±4 


3 


11.5 ±3 


6 


8 


2 


X 


10 11 


±4 


3 


X 


10 10 


1 


2 


X 


10 14 


1 


2 


X 


10 14 


SCSO 


J232342 


3 


-551915 


1 





67 


18.1 ±4 


7 


18.5 ±4 


6 


1 


7 


X 


10 12 


± 1 


5 


X 


10 11 


2 





X 


10 14 


2 


7 


X 


10 14 


SCSO 


J233403 


7- 


-555250 


7 





71 


11.5 ±4 


1 


10.0 ±3 


3 


6 


5 


X 


10 11 


± 1 


5 


X 


10 11 


1 





X 


10 14 


8 


8 


X 


10 13 


SCSO 


J233951 


1 


-551331 


3 





73 


11.6 ±3 


9 


9.9 ±3 


3 


8 


8 


X 


10 11 


± 1 


8 


X 


10 11 


1 


1 


X 


10 14 


1 


3 


X 


10 14 


SCSO 


J233720 


2 


-562115 


1 





75 


10.7 ±4 





7.7 ±2 


9 


5 





X 


10 11 


± 1 


2 


X 


10 11 


8 


5 


X 


10 13 


7 





X 


10 13 



Note. — Catalog of the optical clusters with mass estimates < 3 x 10 14 Mq from the M(Z/2oo) values. Each 
cluster's redshift is the mean photometric redshift computed using the elliptical in the center of the cluster. The ID 
is based on the position of the BCG. 



We have obtained sub-arcsecond astrometric precision 
and sufficient photometric accuracy for the estimation 
of redshift s to Sz < 0.1. Forty- two optical cluster can- 
didates were identified from an area of the sky covering 
^8 deg 2 ; the richness and integrated galaxy luminosity 
of the clusters were measured. Based on correlations be- 
tween these optical observables an d cluster mass as es- 
tablished by SDSS cluster surveys ( Johnston et al.|| 2QQ7 



Reyes et al. 2008) we provide mass estimates (in addi- 
tion to positions and redshifts) for 8 new clusters whose 
inferred masses lie above 3 x 1O 14 M . These clusters are 
all likely to be detected by ACT and SPT if these exper- 
iments reach their expected sensitivity levels. Although 
the uncertainties on the estimated mass and inferred SZE 
signal are large (factors of 2 or so), the accuracy of our 
cluster positions and redshifts are quite good and typical 
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for 4-band imaging survey data. Moreover, it is worth 
noting that the Y — M relation varies quite a bit from 
simulation to simulation (differences of a factor of 2 in M 
at fixed Y"), so our predicted masses are well within the 
uncertainty range of the latest theoretical predictions. In 
an effort to reduce the mass errors we have begun to esti- 
mate weak lensing m asses from the current Blanco data 
QMcInnes et a l. 2009). We also note that significant ad- 
ditional areas of the BCS are now publicly available and 
will be presented in future publications. 

The strength of the SCS is its mult i- wavelength as- 
pect and large sky area coverage. The 23hr region of the 
sky analyzed here has now been surveyed in the UV by 
GALEX and in the X-ray band by XMM-Newton; over 
the next several months as these data become available, 
we will incorporate them into the SCS. Correlation anal- 
yses of these multi-wavelength data should allow us to 
reduce the large errors on inferred cluster masses and 
study the cluster selection biases across wavebands. Se- 
cure confirmation of all candidates, as well as additional 
mass constraints, rest on follow-up optical spectroscopy, 
which we are pursuing at the Southern African Large 
Telescope and elsewhere. Finally, once ACT and SPT 
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